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Abstract 

We show how to compute the coefficients of the double box basis integrals in a massless four- 
point amplitude in terms of tree amplitudes. We show how to choose suitable multidimensional 
contours for performing the required cuts, and derive consistency equations from the requirement 
that integrals of total derivatives vanish. Our formula? for the coefficients can be used either 
analytically or numerically. 



PACS numbers: 11.15.-q, 11.15.Bt, 11.55.Bq, 12.38.-t, 12.38.Bx 



I. INTRODUCTION 



The computation of higher-order corrections to amplitudes in gauge theories is important 
to searches for new physics at modern particle colliders. Next-to-leading order (NLO) correc- 
tions in quantum chromodynamics (QCD), in particular, play an important role in providing 
a reliable quantitative estimate of backgrounds to possible signals of new physics NLO 
corrections to differential cross sections require several ingredients beyond the tree-level am- 
plitudes for the basic process under study: real-emission corrections, with an additional 
emitted gluon, or a gluon splitting into a quark-ant iquark pair; and virtual one-loop correc- 
tions, with a virtual gluon or virtual quark in a closed loop. The required one-loop corrections 
are challenging with traditional Feynman- diagram methods, and become considerably more 
difficult as the number of final-state partons (gluons or quarks) grows. 



The unitarity method [2H18I]. a new method which has emerged over the last decade and 
a half, has rendered such computations tractable. It has made possible a variety of computa- 
tions of one-loop amplitudes, in particular of processes with many partons in the final state. 
In its most recent form, the method can be applied either analytically or purely numeri- 



cally [19H29| . The numerical formalisms underly recent software libraries and programs that 
are being applied to LHC phenomenology. In the current formalism, the one-loop amplitude 
in QCD is written as a sum over a set of basis integrals, with coefficients that are rational 
in external spinors, 

Amplitude = coefficient jlntegraLy + Rational . (1.1) 

j'SBasis 

The integral basis for amplitudes with massless internal lines contains box, triangle, and 
bubble integrals in addition to purely rational terms (dropping all terms of O(e) in the 
dimensional regulator). The coefficients are calculated from products of tree amplitudes, 
typically by performing contour integrals. 

For NLO corrections to some processes, one-loop amplitudes do not suffice. This is the 
case for subprocesses whose leading-order amplitude begins at one loop. An example is the 
gluon fusion to diphoton subprocess, gg — > 77, which is an important background to searches 
for the Higgs boson at the LHC. Although this subprocess is nominally suppressed by a power 
of the strong coupling a s , the large gluon parton density at smaller x can compensate for this 
additional power, giving rise to contributions to cross sections which are comparable to those 
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from tree- level quark-initiated subprocesses [30N32I] . Other examples include production of 



electroweak boson pairs, gg — > Zj, ZZ,W + W~ . NLO corrections to such processes at the 
LHC require the computation of two-loop amplitudes [33 ]. 

Two- loop amplitudes are also required for any studies beyond NLO. Next-to- next-leading 
order (NNLO) fixed-order calculations form the next frontier. The only existing fully- 
exclusive NNLO jet calculations to date are for three-jet production in electron-positron 



annihilation 
LEP 



34j. These are necessary to determine a s to 1% accuracy from jet data at 



351 ]. competitively with other determinations. At the LHC, NNLO calculations will be 



useful for determining an honest theoretical uncertainty estimate on NLO calculations, for 
assessing scale stability in multi-scale processes such as H^+multi-jet production, and will 
also be required for precision measurements of new physics once it is discovered. 

The unitarity method has already been applied to higher-loop amplitudes. At one loop, 
there are different variants of the method. The basic unitarity approach forms a disconti- 
nuity out of the product of two tree amplitudes. Isolating the coefficients of specific basis 
integrals usually still requires performing symbolic algebra on the product of trees; this is not 
well-suited to a numerical approach, and also reduces efficiency of an analytic calculation. 
Basic unitarity corresponds to cutting two propagators in a one-loop amplitude. General- 
ized unitarity cuts more than two propagators at once, isolating fewer integrals. 'Maximal' 
generalized unitarity cuts as many propagators as possible; in combination with contour in- 
tegrals over remaining degrees of freedom, this isolates individual integrals. At higher loops, 
'minimal' generalized unitarity cuts the minimum number of propagators needed to break 
all loops into a product of trees. Each cut is again a product of tree amplitudes, but because 
not all possible propagators are cut, each generalized cut will correspond to several integrals 
and their coefficients, and algebra will again be required to isolate specific integrals and their 
coefficients. This approach does have the advantage of not requiring a basis of integrals. A 



number of calculations have been done t 



lis way, primarily in the M = 4 supersymmetric 



gauge theory and M = 8 supergravity [36H52j|. but including several four-point calcu 



in QCD and supersymmetric theories with less-than-maximal supersymmetry 5314591] . 



at ions 



In this paper, we take the first steps in developing the maximal generalized unitarity 
approach at two loops in a form suitable for both analytic and numerical calculation. We 
show how to extract the coefficient of the planar double box to leading order in the dimen- 
sional regulator e. Higher-loop amplitudes can be written in a similar form to those at one 
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loop (11. ip . as a sum over an integral basis 60], along with possible rational terms. At higher 



loops, however, the coefficients of the basis integrals are no longer functions of the external 
spinors alone, but will depend explicitly on e. Just as at one loop, computing coefficients re- 
quires choosing contours for the unfrozen degrees of freedom. We use the equations relating 
generic tensor integrals to basis or master integrals in order to ensure the consistency and 
completeness of the choice of contours. The extraction of the double-box coefficient bears 
a superficial similarity to the procedure that would be followed in the leading-singularity 



approach 44j, |61(, but unlike the latter, manifestly ensures the consistency of the extraction 



with respect to terms that vanish after integration. Such terms inevitably arise when us- 



ing the integration-by-parts (IBP) approach 62H69j in relating formally-irreducible tensor 



integrals to basis integrals. The extraction of higher-order terms in e or the coefficients 
of integrals with fewer propagators, both of which we leave to future work, would also be 
different. 

During the preparation of this manuscript, a preprint by Mastrolia and Ossola ap- 



peared 79] , analyzing the two-loop integrand in a generalization of the formalism of Ossola, 



Papadopoulos and Pittau (OPP) 11], following a complementary approach to unitarity at 
two loops. 

This paper is organized as follows. In Section [Til we review maximal generalized unitarity 
at one loop, focusing on the computation of the coefficients of the box integral. In Section [TTTl 
we give an outline of the two-loop formalism, and detail the solutions to the cut equations. 
In Section HVl we present the set of constraint equations, and their solutions. In Section [V] 
we give the master formula? for the double-box coefficients, and give some examples of their 
use in Section IVIl We summarize in Section IVIH 

II. MAXIMAL UNITARITY AT ONE LOOP 

We begin by reviewing the derivation of the formula for coefficients of one-loop boxes 



using quadruple cuts, originally written down by Britto, Cachazo, and Feng jg]. We adopt an 
approach and notation that generalize to our derivation for two-loop coefficients in following 
sections. Our starting point is the formal diagrammatic expression for the amplitude, 



i \d Numerator^ (£, • ■ •) • Propagators F (^, ■ ■ ■) 

Ppvnman ■ ' 



(2.1) 



Feynman pi 
diagrams 



where the ellipses represent dependence on external momenta, polarization vectors, and 
spinor strings. Although the whole point of the method is to avoid computing any Feynman 
diagrams explicitly, it is still convenient to refer to them in the abstract, as a means of 
providing the connection to field theory and to Feynman integrals. 

Applying tensor and integral reductions 70|, along with a Gram-determinant identity 
holding through O(e ), we obtain the basic equation (11. ip without any reference to unitarity 
or on-shell conditions. (In a slight abuse of language, we will refer to integrals with no free 
indices, but numerator powers of the loop momentum contracted into external vectors, as 
"tensor integrals".) 

At one loop, it is sufficient for our purposes to concentrate on the four-dimensional com- 
ponents of the loop momentum. (The accompanying integrals must of course be evaluated 
keeping the full (D = 4 — 2e)-dimensional dependence.) In order to derive formula? for the 
coefficients of basis integrals, we apply cuts to both sides of eq. (12.11) . In the basic uni- 
tarity method, we would replace two propagators, separated by a non-null sum of external 
momenta, by delta functions which freeze the loop momenta they carry to their on-shell 
values. In generalized unitarity j^, we would like to apply additional delta functions to 
put additional momenta to on-shell values. However, once we put the momenta carried by 
more than two massless propagators to their on-shell valuess, the solutions to the on-shell 
equations are complex, and taken at face value, the delta functions would actually yield 
zero. 

The same issue arose in the evaluation of the connected prescription 7l| for amplitudes 
in Witten's twistor string theory [72] ; the solution is to use contour integrals instead of 



delta functions 



45 



73 



74j . To do so, we think of complexifying the space in which the 



four- dimensional loop momenta live, from IR 1,3 to C 4 , and taking the integrals on both 
sides of eq. (11. ip to be over a product of contours running along the real axis. We can 
imagine evaluating the loop integrals along other contours as well. New contours that 
will be useful for our purposes are those whose product encircles simultaneous poles in 
all four- dimensional components of the loop momentum. Performing the four- dimensional 
loop-momentum integral over each such contour will yield the residue at the corresponding 
encircled joint or global pole. The residue extracts [] the terms in the integrand which 
contain each of the corresponding propagators, removes the denominators, divides by the 
appropriate Jacobian, and sets the components of the loop momentum to their values at the 
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joint pole. 

The Jacobian is a determinant which arises from the transformation to variables which 
express each denominator factor linearly in a different variable. Unlike a product of delta 
functions, which would produce a factor of the inverse of the absolute value of the Jacobian, 
the transformation here will produce a factor of the inverse of the Jacobian. This ensures 
that the factor is analytic in any variables on which it depends, so that further contour 
integrations can be carried out. 

Notationally, it will still be convenient to use delta functions; to do so, define the product 
of delta functions to yield exactly this contour integral, 

/ t^4i Numerator^, • • -)6(£ 2 )5((£ - h) 2 )5((£ -k x - k 2 ) 2 )5((£ + k A f) = 
J (2vr) 4 

f d A £ Numeratorp(£, ■ ■ •) 

j Tq nt-hyie-h-hyie + hy ' ( • } 

where we have divided out overall factors of 2m associated with each delta function, and 
where Tq is a four-torus encircling the solutions to the simultaneous equations, 

£ 2 = 0, (£-Aa) 2 = 0, {£ - fci - k 2 f = , (£ + fc 4 ) 2 = 0, (2.3) 

and where — in a non-standard bit of notation — we absorb a factor of l/(2iri) into the 
definition of each contour integral, so that evaluating the four-fold contour integral yields a 
sum over residues with no additional factors of 2iri. 

In four-point amplitudes, the external momenta do not suffice to provide a basis for 
arbitrary external vectors; to three of them (say k\, k 2 , and k^), we need to add another 
external vector, for example = e(fi, ki, k 2 , k±). Then we can express all dot products of 
the loop momentum with external vectors in terms of four dot products: £ • k±, £ ■ k 2 , £ • k^, 
and l-v . In reducing integrals, odd powers of v ■£ will give rise to vanishing integrals because 
of parity; even powers can be re-expressed in terms of Gram determinants and thence in 
terms of the other dot products (up to terms involving the (— 2e)-dimensional components 
of the loop momentum). The remaining three dot products can be re-expressed as linear 
combinations of propagator denominators and external invariants, allowing integrals with 
insertions of them in the numerator to be simplified. 

One would be tempted to believe that replacing the original contours running along the 
real axis with some other contour, such as Tq, would leave the equality ( II. II) undisturbed, but 
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this is not quite right, because there are implicitly terms in the integrand of the left-hand side 
that are 'total derivatives', that is terms which integrate to zero. These terms arise during 
the integral reductions described above. They were made explicit in the decomposition of 
the integrand used by OPP . For general contours, the reduction equation will then take 
the form, 

Amplitude = Cjlj + Rat + c j,t u j,t , ( 2 -4) 

j'GBasis j€Basis ^ Total 

Derivatives 

where each Ujj is the integral of an expression Wj t t which would vanish if taken over the 
real slice, for example W\,\ = e(£, k%, k 2 , k$): 

n f d D £ e(£,k u k 2 ,h) 

1 J (2n) D t\t - h) 2 (l -h- k 2 )\i + k A y ' 1 ' ] 

This integral will no longer necessarily vanish if we integrate over general contours in C 4 . 
In this equation, 

e(qi, 92, qa, ^) = e^ a q^qlql , (2.6) 

where s^vXa is the standard antisymmetric Levi-Civita tensor. Another example is the cube 
of the Levi-Civita symbol, = though as it turns out, this latter integrand does 

not contribute any additional equations below. 



FIG. 1: The general quadruple cut in a one- loop box. 

When we perform a quadruple cut, that is the integral over Tq, we will restrict the 
set of Feynman diagrams in the expression for the amplitude to those containing all four 
propagators; cut the propagators; and impose the on-shell conditions corresponding to the 
vanishing of the propagator denominators. If we imagine working in a physical gauge (such 
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as light-cone gauge), this also restricts the cut lines to have physical polarizations. Each 
diagram then falls apart into a product of four diagrams, one corresponding to each corner 
of the box. The sum over diagrams factorizes into a four-fold sum over the tree diagrams 
at each corner, as shown in fig. [TJ Each such sum will give an on-shell tree amplitude, with 
two cut loop momenta and the external legs attached to that corner as arguments, 

d D £ N F 



A i-ioop = J2 [ 

Fey n man jp 



(2n) D D F 



Na N B ^ N c ^ N L 

Fcynmarj 
diagrams 

Ea tree a tree a tree 4tre 
A A A B A C D 



2^ Da 2^ Db 2^ Dc 2^ Dd y z - l > 

hclicitics Fcynman Fcynman ^ Fcynman q Fcynman j-y 



species diagrams diagrams diagrams diagrams 

i tree 



hclicitics 
species 



(If we had not initially used a physical gauge, it is at this stage, summing over all diagrams, 
that we would recover the restriction to physical polarizations.) Integrating over Tq in 
eq. ( 12. 4p will give us the following equation, 

53 ATA^A^A^ = c box + ^ c' hoXtt U hox , t (2.8) 

hclicitics £ 
species 

where the Jacobian has canceled out of the equality, and the sum on the right-hand side 
runs over possible total derivatives with a box integrand. In order to solve for the desired 
coefficient Cb ox , or equivalently to ensure that the equality in eq. (11. ip is maintained, we 
must evaluate the integral over a linear combination of new contours such that all possible 
integrals of 'total derivatives' Ut are projected out. As we will show later, this requirement 
gives us constraints that determine the allowed combinations of contours, and in turn the 
equations for the coefficients of the box integrals. 

In the case of the one-loop box integral, the joint-pole equations are given by eq. (12 .3p or 
equivalently by, 

£ 2 = 0, £-ki = 0, t-k 2 = S -, £-k 4 = 0, (2.9) 

which form makes it clear that there are two distinct solutions, and hence two distinct 
contours. The one 'total derivative' we must consider is the e expression U\ given above in 
eq. (12. 5p . It turns out that it evaluates to compensating values on the two solutions, so that 
summing over them projects it out, and hence gives an equation for the coefficient of the 
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box in terms of the product of tree amplitudes at each corner, 

c b ox = ^ AT AT AT AT- (2.10) 

hclicitics 
species 

As an example, study the coefficient of the one-mass box with m\ 7^ 0. Parametrize the 
four-dimensional part of the loop momentum as follows, 

e = a lK + a# + J 0L_ (I", r |2-> + (2-| * . (2.11) 

A general contour integral for the four- dimensional part of the box integral then takes the 
form, 

<i 4 a,- 1 



c J a {a.10,2 — uia^a4){aia2 — uia^a^ — a 2 )(a\a 2 — uja^a^ — a 2 — «i + 1) 

x- ] (2.12) 

(«i«2 — uia^a^ + a\tj s + a 2 u/s + a 3 + a i) 

where u = s 2 /(tu). In this expression, J a is the Jacobian that arises from changing variables 
from the to the aij. (We do not need its explicit form, only the knowledge that it is 
independent of the a^, a consequence of the linearity of the in the a«.) 
The cut equations ( 12.91) then take the form, 

(aia 2 - ^-a 3 a^js = 0, 

a 2 s = , 

ais = s, (2-13) 

a\t + a 2 u + + a^s = , 

which have two solutions, 

51 : a\ = 1 , a 2 = , «3 = — , «4 = ; 

s 

5 2 : «i = 1 , a 2 = , a 3 = , a 4 = — . (2-14) 

s 

If we define Cj(v) to be a small circle in the complex <x,-plane that encloses the point v, 
then the two contours we must consider are, 

Tt = d(l) x C 2 (0) x C 3 (-t/s) x C 4 (0) ; 
and T 2 = d(l) x C 2 (0) x C 3 (0) x C 4 (-t/s) . (2.15) 



We can evaluate the four- fold integral (12.121) by 'global residues' |45|, |75|]. The sign of 
the result will depend on the orientation chosen for the contour; but this sign will drop out 
of final formulae for integral coefficients so long as this orientation is chosen consistently 
throughout the calculation. To do so, we should first change to variables where each pole 
is in a different variable, and where the denominators are linear in that variable with unit 
coefficient. The Jacobian from this change of variables will take the form, 



where 



i,j \OOii 



fl = O1O2 — UJOtzOLi , 

f 2 = Oi0 2 - ua 3 a4 - a 2 , 

fs = a\a 2 — uia^a^ — a 2 — oi + 1 , 

f± = atict 2 — WO3O4 + (X\tj s + a 2 u/s + 0:3 + 0:4 



(2.16) 



(2.17) 
(2.18) 



We find, 

Ji = w(o 4 -o 3 ). (2.19) 

Evaluating the box integral with a numerator Num(£, • • •) along a contour given by a linear 
combination of the two with weights <2j, we obtain, 



s-^J' 1 ( at JT^Num^, • • •) + a 2 Jr 1 Num(£, • • •) ] 
V <Sj s 2 J 



(2.20) 



Using the parametrization (12.111) . we find the following expression for the Levi-Civita symbol 
we need, 

e(£, fa, k 2 ,k 4 ) = a 3 se f -^j^p^ , fa, k 2 , fc 4 ) + a * s£ ( %0$[T\ » k ^ k *) 

= s(a 3 -ai)e ^l { }l^ ,fa,k 2 ,kA . (2.21) 

The constraint that U\ = on the quadruple cut then implies that, 

-s-*J- 1 cj- 1 (a 1 + a 2 ) = 0. (2.22) 

so that a 2 = —a\. Higher odd powers of the Levi-Civita tensor lead to the same constraint. 
If we evaluate both sides of eq. (11.11) on the linear combination of contours, we find, 
2 4 



tree 



i=l 



s 4 J a 1 c I OiJf 



Si 



r-1 



Si 



a 2 J[ 



-1 



s 2 



(2.23) 



10 



where the product is over the tree amplitudes associated to each of the four vertices of the 
quadruple-cut box integral in fig. [TJ Substituting in the solution to eq. (12. 22 p . we find for 
the coefficient of the one-loop box, 

u 2 4 

C-box 



2 4 

Hll A T e , (2-24) 



2 

i=i i=i 



which is just eq. f )2.10p when summed over possible helicity assignments and species circu- 
lating in the loop. 

In the following sections, we show how to generalize these considerations to two loops. 



III. MAXIMAL CUTS AT TWO LOOPS 



Our basic approach to the planar double box at two loops will be similar to that reviewed 
above at one loop. We use a convenient parametrization of the loop momenta and choose 
new contours of integration to freeze the momenta flowing through all propagators. We 
choose those contours so that constraint equations arising from consistency conditions are 
satisfied. Unlike the procedure at one loop, cutting all seven propagators does not freeze 
all components of both loop momenta, so we must choose new contours for the remaining 
unfrozen degrees of freedom as well. In addition, we have a much larger and richer set 
of consistency conditions arising from IBP identities. Once we have solved the constraint 
equations, we will solve for the coefficients of specific basis integrals. 

k 2 fa 




k\ fa 



FIG. 2: The double box integral P 2 ** 2 . 

In this section, we give a convenient parametrization of the loop momenta, and use it to 
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solve on-shell equations. We list these solutions below, along with the poles and possible 
contours for the remaining unfrozen degrees of freedom. 
The two-loop double box integral, shown in fig. [21 is, 

d D e 1 d D i 2 1 



R 



2,2 



(2tt)d (2tt) d £1(1, - hYih - K X2 )\£ X + £ 2 )H 2 (£ 2 - K)\£ 2 - K u ) 



(3.1) 



where Ki...j = k t + ■ ■ ■ + kj, and the notation follows ref. [60]. 



We will focus in this paper on extracting coefficients of basis integrals only to leading 
order in the dimensional regulator e, for which it suffices to consider the four- dimensional 
components of the loop momentum as far as cuts are concerned. The double box has seven 
propagators; if we cut all of them, that is put all of the momenta they are carrying to on-shell 
values, we will be left with one additional degree of freedom. To cut the momenta in this 
way, we must shift the contours of integration for the components of the two loop momenta 
l\ and £ 2 to encircle the joint solutions to the on-shell equations, 







(£i-h) 2 = 0, {^-K^f 







(£ 2 -^ 3 4) 2 = 0,(£ 1 +£ 2 ) i 







£ 2 2 = 0, (£ 2 -hf 







(3.2) 




FIG. 3: The heptacut double box. 

As explained in Section [Til we can write the four- dimensional heptacut integral symboli- 
cally as, 



d 4 £i d% 2 



l)6((£ 1 -k i y)5((£ 1 -K 12 y)5((£ 1 + £ 2 



(2tt) 4 (2tt) 4 

x5(£ 2 2 )5((£ 2 - h) 2 )5((£ 2 - K u ) 2 ) , (3.3) 
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again dropping overall factors of 2iri associated with the delta functions. This heptacut is 
depicted in fig. [3j 

To solve the on-shell equations, we use the following parametrization of the loop momenta, 

Si2a3 uln _> s 12 a A 



a±ki + a 2 k^ + 



2 (1 4) [4 2] 
S12A5 



<3~| 7" |4~) + 



2 (2 4) [4 1] 

Sl 2 /?4 



(2-irii-> 

(4-| 7 M |3-) 



(3.4) 



2(41) [13] 

Using this parametrization, the six corresponding heptacut equations involving only one 
loop momentum are, 



{£1 - hf 



(£i-K 



12 j 



(£2 - K M f 



S\2 
S\2 
S\2 

S12 
S12 



ctia 2 + 



x(x+i) 

O; - 1)« 2 + 











x(x+i) 

(«i-l)(«2-l) + ^y 
" *) + jgfe) = 

(A-i)CA-i) + 5g&j 



(3.5) 



where 



X = 



£14 
S12 



(3.6) 



We can simplify these equations, obtaining 



«1 = 1 , «2 = , «3«4 = 

/?! = (), & = /3 3 /3 4 = . 



(3.7) 



These equations have four distinct solutions. If we substitute these values into eq. (I3.2p . we 
find for the last equation, 







(£ l +£ 2 f = 2£ l -£ 2 



2 fc? + 



2(14) [4 2] 



(l-| 7 l2-) + 



Si2«4 



2 (2 4) [4 1] 



(21/in 



X, ^ + 2(3T)[T4] ( hJ ) + 2741)113] ( l7j } 



(3.8) 



For two of the four solutions to eqs. ( 13. 7p . this equation has two solutions, so that overall 
we find six solutions to the heptacut equations (13.51 13.8j) . To each of the six solutions Sj, 
we can associate a seven-torus in the parameters ctj and fa that encircles the solution. 
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For the solution a 4 = = /3 4 , the last equation (13. 8 p simplifies to, 



°=H+ it?) ( <i4) -w)' (a9) 

which has two distinct solutions, 

S 1 : « 3 = -x , /3 3 arbitrary ; 

<$2 : = — X ) a 3 arbitrary. (3.10) 

In all solutions, we will change variables so that the remaining degree of freedom is called z. 
Likewise, the solution a 3 = = /3 3 also yields two solutions to eq. (I3.8p . 

5 3 : a A = -x , P>4 = z\ 

54 : 04 = — x , 0*4 = z . (3-11) 

For the remaining two solutions, the last equation (13. 8 p does not factorize, and we obtain 
only one solution; for a 3 = = /3 4 , 

S 5 : a4 = z } < 9 3 = -( x + l) Z + X ; (3.12) 

z + X + 1 



and for a 4 = = /3 3 



5 6 : « 3 = z, /3 4 = -( X + l) ■ (3.13) 

z + X + 1 



In the last two solutions, we could equally well have chosen a different parametrization, 
where (3 3 or /3 4 respectively are set to z. This just amounts to a change of variables, of 
course, but does break the manifest symmetry between the two loops. 

The existence of six kinematic solutions can also be understood from holomorphicity 
considerations of the cuts. When we cut all propagators, each of the six vertices in the 
double box has three massless momenta attached. We can write these momenta in terms 
of spinors, = X a a l ^ (k X a . Momentum conservation at each vertex 7^ then implies that 
either, 

1) the holomorphic spinors A of the momenta are collinear (proportional), A a oc A?, oc A c . 
We will depict such a vertex using a circled plus (©). Such a vertex would allow only 



an MHV tree amplitude to be attached (of course the holomorphicity properties of the 
cut are independent of any tree amplitude). 
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Solution Si, obtained by setting 

«3 = -X , fa = z, 




Solution 1S2, obtained by setting 



«4 = , 



ft = 0. 



«3 = Z , 
«4 = , 



ft = -X : 
04 = 0. 




Solution 53, obtained by setting 

«3 = , ft = , 

04 = -X , ft = 2 . 




Solution 54, obtained by setting 

as = , ft = , 
ai = z, ft = -x ■ 




Solution 55, obtained by setting 
a 3 = 0, ft = -( x + l)_!±X_, 




«4 = Z . 



ft = 0. 



Solution Sq, obtained by setting 

«3 = Z , ft = , 

a 4 = 0, ft = -(x + l)ifgi- 



FIG. 4: The six solutions to the heptacut equations for the two-loop planar double box. 
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2) the antiholomorphic spinors A of the momenta are collinear, A a oc A;, oc A c . We will 
depict such a vertex using a circled minus (0). Such a vertex would allow only an 
MHV tree amplitude to be attached. 

For general kinematics, neither the external holomorphic spinors Aj nor the external anti- 
holomorphic spinors Xj are collinear. A configuration with an uninterrupted chain of either 
© or vertices connecting any two external legs is thus disallowed. There are exactly six 
ways of assigning these two labelings to vertices avoiding such chains, hence six solutions. 
The six solutions are shown diagrammatically in fig. HI (The labeling of holomorphically- 
collinear vertices as 0, and of antiholomorphically-collinear ones as is not uniform in the 
literature.) 

In evaluating the contour integrals represented by the delta functions in eq. (13. 3p . we 
encounter two Jacobians: one from changing variables from the components of Ej to the 
and Pi] and one from actually performing the contour integrals in the latter variables. It 
is the latter Jacobian that is important for our purposes. The former Jacobian is equal to 
J a J/3, where 

J Q = detffi= - f 12 - , Jg = detgf = - f 12 - ■ (3.14) 

To evaluate the latter Jacobian, we may note that three of the delta functions (or equiv- 
alently three of the contour integrals) involve only a variables, and three involve only (3 
variables. We can thus split up the problem into three steps: computing the Jacobian asso- 
ciated with £i, that is with the a variables alone; computing the Jacobian associated with 
£2, that is with the variables alone; and finally, computing the Jacobian associated with 
the middle propagator, involving both £1 and £ 2 - 

For each of the six solutions, we must compute the Jacobian independently. As an 
example, consider the second solution <S 2 - The first Jacobian arises from considering the 
integral, 



J doL\_doLidoL^b s i2 (0:10:2 + x °x+i) ) ^ s i2 ((ai - l)a 2 + 



Q3 "4 



x8 



Su (( ai -l)(a 2 -l) + 3g^)] , (3.15) 
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associated with the i\ loop. Define 

( 9i(a^ 

92(012 
^3(04) J 
The first Jacobian is then, 



( s 12 ( ai a 2 + 



S12 



(3.16) 



/ 



Ji = det 

m aa,- 



sf 2 det 



Q3 \ 



a 2 



Ql 



«2 «1 



x(x+i) 
x(x+i) 



'12 



(3.17) 



(As explained in Section [Til the Jacobians will appear in the denominator as determinants 
rather than as absolute values of determinants.) Similarly, the second Jacobian arises from 
considering the integral, 



x5 



S12 /5i(/5 2 -l) + 



x(x+i) 



associated with the £2 loop. Define 
^i(/3i^ 

\h 3 ((3 4 ) J 

The second Jacobian is then 



si 2 ((/3i-l)(/3 2 -l) + ^) 



sis f/3i(/3 2 -l) + ^ 



(3.18) 



x(x+i) 



(3.19) 



/ 



J? = det 



J 30 j 



s 12 det 



A 0i 

02-1 Pi 



V 



ft \ 

x(x+i) 

§3 

x(x+i) 

ft 



„3 
6 12 



(3.20) 



^ 2 - 1 & - 1 x(x+i) / 
The remaining integration we must consider is over a.3 and 0s 

S12 



2 j Ji-h 

daisdfiz 



32s 2 2 



2X 
5 



(a 3 + x)(/3 3 + x) 



^(«3 + x)(/?3 + x) 



(3.21) 



which leaves a remaining contour integration over z {i.e. a^), along with the overall inverse 
Jacobian, 

J~Hz) = — 1 . (3.22) 



16s? 2 z(z + x) 
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The computation for the other five solutions is similar; it turns out that we obtain the same 
overall Jacobian for all solutions. The contour for the z integration remains to be chosen; 
for this solution, there are two possible non-trivial contours, one encircling z — 0, and the 
other, encircling z = —\. (We set aside a possible non-trivial contour encircling z = oo, 
as its contribution when integrating an arbitrary multiplying function f(z) sums to zero 
when combined with the contributions of these two contours.) The pole at z — —\ is the 



eighth pole in the octacut of ref. 6l|- In addition, for solutions £5^, the denominator of 
(eqs. (13.121) and (13.131) ) can give rise to additional poles at z = ~x — 1 in tensor integrals. 
(As noted in Section [III in a slight abuse of language, we refer to integrals with no free 
indices, but numerator powers of the loop momenta contracted into external vectors, as 
"tensor integrals".) 

Collecting the information above, we have the following contours we can utilize in seeking 
equations for integral coefficients, 



Ti,i 

Tl,2 
T2,l 

T22 
^3,2 

^4,2 

Tw 
T 5 , 3 



T xC a ,(- 
T xC a3 (- 

To X C a3=z 

To x C as=z 
T x C Q3 (0 
T xC Q3 (0 
T xC Q3 (0 
T x C Q3 (0 
T x C Q3 (0 

T x C Q3 (0 

T x C Q3 (0 

To x C as=z 

To x C as=z 

Tq x C as=z 



x)xC fl4 (0)xC ft=2 (0)xC ft (0), 
x)xCJO)xC ft=z (-x)x^(0), 
0)xC Q4 (0)xC fe (- X )xC /34 (0), 
-x)xQ0)xC ft (- x )xC ft (0), 
xC a4 (0)xC ft (0)xC ft=2 (0), 

xC a4=z (0)xC ft (0)xC^ 4 (-x), 
x C a4=z (-x) x Cp 3 (0) x C A (-x) , 
x C ai=z (0) x Cp 3 + *)(* + *) ) x c- (o) 



(3.23) 



x C ai=z {-x) x C Ps 



z + X + ± 
z + X+l 



x c a4=z (- x - 1) x c p , A [ - {1 t x){z+ , x) ) x c P M , 



0)xQ)xC ft (0)xC ft 
-X) x C Q4 (0) x C fe (0) x C P4 
-X-l)xCJ)xC fe (0)xC ft 



z + X + l 
(1 + X)(z + X) 
z + X + ± 
(1 + X)(z + X ) 



z + X + ^ 
{1 + X){z + X) 
z + X + 1 
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where each subscript denotes the variable in whose plane the circle lies, and where 

T = C ai (l) x C Q2 (0) x C h (0) x C fe (l) , (3.24) 

corresponding to the on-shell values in eq. (13. Tp . We will call the complete contours, including 
a choice of contour for z, the 'augmented heptacut'. 

Naively, we could deform the original contour of integration for the double box ( 13. ip . 
along the product of real axes for all components of l\ and £ 2 , to any linear combination 
of contours in eq. (I3.23P that we wish. However, an arbitrary deformation will not preserve 
the vanishing of total derivatives, analogs to U\ given in eq. ( 12. 5p . In order to ensure that 
such objects vanish as they must, we impose constraints on the contours. We derive these 
in the next section. 



IV. CONSTRAINT EQUATIONS FOR CONTOURS 

Integral reductions are implicitly part of the simplifications applied to a sum over Feyn- 
man diagrams in order to obtain the basic equation at either one loop (II. ip . or at two 
loops, 

Amplitude = coefficient j Integral^ + Rational . (4-1) 

j'£Two-Loop Basis 

The basis at two loops will contain integrals with up to eight propagators in the planar 
case [60], though a specific complete and independent choice of integrals for a general am- 
plitude has not yet been written down. (The same restriction to eight propagators or fewer 
presumably applies in the non-planar case as well, using arguments along the same lines as 
given in ref. [60].) 

As we saw in Section m integral reductions at one loop involve only rewriting dot products 
of the loop momentum in terms of linear combinations of propagators and external invariants, 
along with the use of Lorentz invariance and parity to eliminate some integrals. For the box 
integral, in particular, the only non-trivial constraint arises from the use of parity, which 
requires that 

d D £ e(£,k u k 2 ,h) 



(2ti) d £ 2 (£- k x f{£- kt - k 2 y{£ + h) 2 °' ^ 
This constraint must be respected by the unitarity procedure; otherwise, applying a cut 
to the original integral and to the integral after reduction would yield different, and hence 
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inconsistent, answers. At one loop, it gives rise to one constraint equation, which fixes the 
relative normalization of the contours encircling the two solutions to the on-shell equations. 

Similar constraints arise at two loops, though we have a greater variety of Levi-Civita 
symbols to consider. Denoting the insertion of the function f(£i, £ 2 ) i n the numerator of the 
double box by 

d D £ x d D £ 2 f(£ h £ 2 ) 



(2ir) D (2n) D £\{£ x - k x f{£ x - K 12 )\£, + £ 2 )H\{£ 2 - k,)\£ 2 - K u f ' 

(4.3) 



we must require that the vanishing of the following integrals, 

P^Hh, k 2 , k 3 , h)) , p;:M£ 2 , k 2 , k 3 , h)) , P%Wuh, h, k 2 )} , 

P**M£i,£2,h,k 3 )}, and P;; 2 [e(^^,h,k 3 )}, (4.4) 

continues to hold for integration over our chosen linear combination of contours. This is 
the complete set of Levi-Civita symbols that arises during integral reduction, after using 
momentum conservation. 

At two loops, additional reductions are required in order to arrive at a linearly- 



independent set of basis integrals. These are usually obtained through IBP relations 62 



63 



66M68|; that is, they correspond to adding a total derivative to the original integrand. Each 
such total derivative, or equivalently each non-trivial reduction identity, gives rise to a con- 
straint requiring that the unitarity procedure give vanishing coefficients for the additional 
terms; or equivalently, that the unitarity procedure respect the reduction equations. This 
is not automatically true contour-by-contour, and hence gives rise to non-trivial constraints 
on the choice of contours, and the weighting of different solutions. 

In two-loop four-point amplitudes, we can express all dot products of loop momenta 
with external vectors in terms of eight dot products: £j ■ fci, £j ■ k 2 , £j ■ k±, and £j ■ v, where 

= e(/i, k\, k 2 , ki). Just as at one loop, odd powers of v will give rise to vanishing integrals, 
as expressed in the Levi-Civita constraints discussed above. Even powers can again be re- 
expressed in terms of the other dot products (up to terms involving the (— 2e)-dimensional 
components of the loop momentum). All integrals can then be rewritten in terms of the six 
dot products of the loop momenta with the external momenta. 

Of these six dot products, three of them — £ 1 ■ ki, £\ ■ k 2 , £ 2 ■ fc 4 — can be rewritten as 
linear combinations of the propagator denominators and external invariants. One additional 
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dot product of £2 — say £2 • ^2 — can be rewritten in terms of the remaining two (£1 ■ and 
£2 ■ hi), propagator denominators, and external invariants. The remaining two dot products 
are called irreducible. At a first stage, then, before using IBP identities, we can reduce an 
arbitrary double-box integral appearing in a gauge-theory amplitude to a linear combination 
of the 22 different integrals that can arise with powers of the two irreducible numerators. 
We have the following naively-irreducible integrals, 

p 2 ** 2 [i],p 2 *; 2 [^-4], p 2 *; 2 pi-4) 2 ], A*,*2pi-4) 3 ], p^ih-ti], 
p%[{h ■ £ 2 ) (h ■ £1% p;; 2 [(h ■ £2? (h ■ £1)], p;; 2 [(h ■ e 2 ) 3 (h ■ h)l 

P 2 %(h-£i) 2 }, P^Kh-Q P 2 * 2 pi • £ 2 ) 2 (h-^n 

P&iki ■ 4) 3 (h ■ £1)% P;; 2 [(h ■ £1)% P^iih ■ 40 (h ■ £i)% (4.5) 
P%[(h ■ £ 2 ) 2 (h ■ £,)% P^[(h ■ £ 2 f (h ■ £,)% P^Kh ■ £i)% 
p%[(h ■ £ 2 ) (h ■ £1)% P^Kh ■ £2? (h ■ £1)% P%[{h ■ £2)% 
P^Kh ■ £ 2 f (h ■ £1)], P 2 %(h ■ £ 2 ) 4 (h ■ £ x f] . 

In the massless case, it turns out that there are 20 IBP relations between these integrals, 
which allow further reductions. These reductions allow us to pick certain pairs, for example, 

P** 2 [l] and P;* 2 [li-h], (4.6) 

as master integrals for the set in eq. (14. 5p . and thus also as basis integrals for an amplitude. 
The remaining integrals are given in terms of these two by linear equations, for example 

p;; 2 [£ 2 -h} 

p^Kh-h^-h)} 
P&fr ■ h) 2 } 

P2W1 ■ k 4 )(£ 2 ■ h) 2 ] 
P 2 * 2 [(4 • k A ) 2 {£ 2 ■ h)] 



--p;* 2 [£ 1 -h] + 
2.: 

i . (x + 3e)si2 



-\xs\ 2 P^[l]-\s l2 P^[h-h] + 



4(1 - 2e) Z ' 2L 2(1 - 2e) 



P2I [i] + , ' v p 2*l[h ■ h] + 



4(1 - 2e) ^ L 1 2(1 - 2e) 
< ' - 3e)xs? 2pr2[1] + (3-9e-2e X )s 2 2 p ^ [h ^ + ^ 



16(1 -2e) ^ L J 8(1 -2e) 

"^2,2 [1] + ^ 7T\ P 2, 2 R • k 4 + 



16(1 -2e) J 8(1 -2e) 

6x(l -2x- 3e)sf 2 , 
16(l-e)(l-2e) 2 ' 2[ 1 
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(2 x 2 -3e(l-2 X ) + e 2 (9 + 2 X )) g 2 
+ 8(l-e)(l-2e) 2 < ai 1 ' 4 ' 



P**M-k 4 ) 3 ]= ex{ ] 2 * 3e) 1 2 P 2 ** 2 [l] 
2,2LV 4; j i6(l- e )(l-2e) 2,21 J 

(2y 2 - 3e(l - 2y) + e 2 (9 + 2y))s 2 „^ rjl , . 

+ v_t ^ a; t v ~ A;; 12 -P 2 *o • A; 4 + • • • , 

8(l-e)(l-2e) 2 ' 2L J 

where the ellipses denote additional integrals with fewer propagators. We must require 
that these equations (and the other 12 we do not display explicitly) are preserved by the 
choice of contours. The contour integrals which implement the augmented heptacut will 
yield vanishing results for the integrals with fewer propagators, so they do not enter the 
constraint equations. As we are considering only four-dimensional cuts, the augmented 
heptacuts are effectively four-dimensional. 

In order to find the explicit form of the constraint equations, denote the weight of contour 

Pj,r by dr,j i 

aij — > encircling z = for solution Sj , 

02j — > encircling z = — x for solution Sj , (4.8) 

a 3 j — > encircling z = —\ — 1 for solution Sj . 

For a numerator insertion of f(£i, £2) in the numerator of the double box, the augmented 
heptacut is then, 

4 2 



par am 



/] / J ar >3 f d 4 ®id 4 f3i fihAi) x Propagators(^i, £2) 

j=l r =l J T j, r 
6 3 „ 

+ VVa rj i d*aid% f(£ij 2 ) x Propagators (4 ,4) , (4.9) 

' ' — J It param 

j=5 r =l J 1 o,r F 

where the notation | param indicates that we use the parametrization of £\ and £2 given in 
eq. ( 13 .41) . The signs in front of each coefficient a r j in the result will depend on the orientation 
chosen for the corresponding contour; but this sign will drop out of final formulae for integral 
coefficients so long as this orientation is chosen consistently throughout the calculation. 

We can write down a compact expression for the augmented heptacut of the general 
tensor integral, 

1 



128 
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/ Sl2 Nm-3 f d,Z 



/ Sl2 \m+n-3 /• (fz 



+ -fa , 

V 2 / / r6 z V^ + X + l 



(4.10) 



where Tj denotes the z component of ^ r a r jTj >r , and where in our notation, the contour 
integral implicitly includes a factor of l/(27ri), as noted in Section HT1 

We can evaluate this expression using the contours as weighted in eq. (14. 9ft ; we find, 

1 6 

p 2 *; 2 [(4 • h) m ] = (^) m_1 £ fl u ■ ( 4 - 12 ) 



'12 



3^1,3 



P 2 ** 2 [(£ 2 • fci) n ] = ^-L- (^) n 1 (-02,6 + a 3 , 6 - a 2j5 + a 3l5 + a lfl + a 1>3 ) , (4.13) 

1 / a, „ \ m+n— 2 r 

P*; 2 [(l! ■ h) m (£ 2 ■ hT] = -— (-- [0 (min(m, n) -§)(* + 1)(* + 2) 

+ (min(m, n) — |) (m + n — 3)(x + 1) + G (min(m, n) — |) 

x(a 3 ,6 + a3, 5 ), (4.14) 

where m, n > 1 and the last result is valid only for < m + n < 6 and < m, n < 4 (corre- 
sponding to the numerator insertions allowed in gauge theory in D = 4 — 2e dimensions). 

With these expressions, we now turn to the constraint equations. Let us begin with the 
equations arising from the insertion of Levi-Civita tensors (14. 4p . Start with £(£%, k 2 ,k3, k^j, 

o = P^WuhMM)] 

, e ( S 12 (IZhUPll ^ b b b\ , c f b^ 1 -L Sl2 (1~|7 M |2~) t, 7 , 

() ... cfo 6 ^! 2 (l-|4|2-> X,»2,K3,K4J ^ £ dz £ \ K ^ + 2 (F |4|2~> Z ' K 2, K 3, 



16s? 2 \ / ri z z + x ' Jt 2 z z + X 

r cfb^ S 12 (2~l7 M |l~) b b b\ , c ( b^ 1 _l_ 512 ( 2 ~|7 M |1~) y jL 7, J, \ 

2~ <2-|4|l-) X, K 2,K3,K±J [ dz £ + ~ 1?=W^Y ' 2 ' 1 A ) 

JTz z z + X Jr 4 z z + x 

, p-fb^-L^ (2~|7 M |1~) b b b\ „ , e ( b^ -L S 12 (1~|7 M |2~) j. t, j. \ 

2 (2-141!-) Z,K 2 ,K3,K4j fdz £ \ K l + 2 <1" |4|2"> Z ' K 2, «3, ^4 j 

r 5 z Z + X Jr 6 z z + x 
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(4.15) 



Evaluating this expression on the augmented heptacut (14.91) . we obtain, 
1 



32s 2 2 



f{l~\^\2~) 

(a 2) 2 + «2,6 — a l,l + a 2,l) £ ( ^ |2_^ ) ^3) ^4 

/(2"| 7^ II") 

+ (a 2 ,5 + «2,4 - «1,3 + «2,3) £ ( /r ,_| TTT^N > k 2 , &3, &4 



(2"|4|l-> 

1 /(1~|7'*|2-) 

- ^"(02,2 + «2,6 — a 2,5 — a 2,4 _ 0-1,1 + 2j l + a l,3 — ^2,3)^ [ ,-,_,„ in-\ > ^2, &3> &4 



(4.16) 



•'W.s-f/ - ' - -' ~* ' ^ ' - V (l-|4|2-> 

where the last line follows from the fact that the two Levi-Civita symbols appearing on the 
first line are equal but opposite in value. 

Similarly, from the insertion of e(£2, k 2 , &3, k^) one finds 



= P** 2 [e(£ 2 ,k 2 ,k 3 ,h)} 

1 . . 

- a\ t 2 + 0,2,2 + 0>1,6 — «3,6 — Ol,5 + a 3,5 + Ol,4 — a 2,4 + «2,1 — a 2,Z) 







32s 2 
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(3"|7 M |4~) 

X£| - , , , , ; , fc 2 ,fc 3 ,fc 4 



(3-|l|4-) 

and from the insertion of s(£i,£ 2 ,ki,kj) with (2, j) e {(1,2), (1,3), (2,3)} one finds 



= P;* 2 [e(£ 1 ,£2,k 1 ,k 2 )} 
1 







32s 2 2 



(o 2 ,6 ~~ a 3,6 — a 2,5 + a 3,5 — a l,l + a l,3) 8 



(1-1 




2-) 


(1- 


|4|2-) 



^2; ^3) ^4 



o = P 2 ** 2 [£(4, £2,h,k 3 )} 
l 



o 



(l"|7^|2~) 

2 ( a 2,6 — a 2,5 — a l,l + ai,3)£ ( /-,_!, ) &2, ^3; ^4 



32 s -^2 

1 







32s 2 2 



(l-|4|2-> 

(oi, 2 — 01,6 + a 2,6 + 0^5 — a 2 ,5 — ax,4 — 0^1 + 0^3) 



(1~|7 M 12") 
xel , , , , k 2 , fc 3 , fc 4 



(1-|4|2"> ' 

These results combine to give the constraints, 

G 2,2 + O2.6 — 02,5 — 2j 4 — Oil + 02,1 + Oi ; 3 — <2 2 ,3 = 

Ol )2 — 0-2,2 ~~ Ol,6 + 03,6 + a l,5 — 03,5 — Ol,4 + °2,4 — 2j i + 0, 2j 3 = 

0-2,6 ~~ a 3,6 — a 2,5 + 03,5 — a l,l + a l,3 — 

02,6 — °2,5 — Ol,l + 01,3 = 

«1,2 - «1,6 + 02,6 + «1,5 - «2,5 - «1,4 ~ «1,1 + Ol,3 = 
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(4.17) 



(4.18) 



(4.19) 



or equivalently, 



01,2 



Ol,6 + °1,5 — Ol,4 



^2,2 



«2,4 + ^2,1 — a 2,3 



a 2,6 



a 2,5 — °1,1 + a l,3 







(4.20) 



^3,6 — a 3,5 



0. 



This set has one equation less: not all the equations from the Levi-Civita symbols are 
independent. We see that these equations are solved by insisting that the contours for 
complex-conjugate pairs of solutions (Si 4-> S 3 , S 2 O S±, and S 5 -H- S 6 ) carry equal weights. 
This nicely generalizes the one-loop constraint on contours. However, these are not the only 
possible solutions; solutions which do not insist complex-conjugate pairs carry equal weight 
are also possible. 

We next impose the constraints following from the IBP reductions. Evaluating both sides 
of equations (14. 7ft along with the remaining 12 reduction equations not displayed above, and 
setting e = 0, we find two additional constraint equations, 



In principle, one might expect 18 additional equations from the remaining reduction iden- 
tities; but these all turn out to be automatically satisfied on the solutions of this pair of 
equations. 

Beyond ensuring that all the reduction identities are valid, we ultimately want to deter- 
mine the coefficients of the two basis integrals (14.61) . Because the system of equations leaves 
many undetermined weights a r j, we have the freedom to choose values which also kill one 
or the other of the basis integrals. That is, we can choose contours for which one or the 
other of the basis integrals has vanishing augmented heptacut. To project out the second 
basis integral, -P 22 [^i ' we should also require that eq. (I4.12p vanish, 



a%,2 + a 1>6 + ai j5 + a M 



-a 2 ,6 + Ct3,6 - Q-2,5 + 0-3,5 + 0-1,1 + 0,1,3 , 




(4.21) 




(4.22) 



m,3 



To project out the first basis integral, -P 2 * 2 [l], we should require that eq. (14.1 ip vanish, 



6 




(4.23) 
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The following values, 



ai,i = 


-2u + v , 


02,1 


= u, 




a>i,2 = 


-2u + v, 


«2,2 


= u, 




«1,3 = 


-2u + v, 


02,3 


= U, 




Oi 5 4 = 


-2u + v, 


02,4 


= U, 




ai,5 = 


2u — v , 


02,5 


= v, 


a 3j5 = 2u, 


ai,6 = 


2u — v , 


02,6 


= V, 


a 3j6 = 2u, 



(4.24) 



(where u,v are real parameters) solve all the constraint equations ()4.20[ 14.21)) . and also set 
the basis integral P 2 *2^i ' ^4] ^° zero, thereby allowing us to extract the coefficient of the 
first basis integral, f^W- We will call a specific choice of contours weighted by these values 
Pi, leaving the dependence on u and v implicit. A particularly simple solution is given by 
u = \ and v = 1. 

Similarly, the following values, 



Ol,l = 


-2u + v, 


02,1 


= u, 






Ol,2 = 


—2u + v , 


02,2 


= u, 






Ol,3 = 


-2u + v, 


«2,3 


= u, 






Ol,4 = 


—2u + v , 


02,4 


= u, 






0-1,5 = 


6u — v , 


02,5 


= V, 


03,5 


= 6u , 


a l,6 — 


6u — v , 


02,6 


= 


03,6 


= 6u , 



(4.25) 



(where again u, v are real parameters) solve all the constraint equations (I4.20[ POT]) , sets to 
zero the basis integral i^fl]; an d thereby extracts the coefficient of P 2 *2^i ' ^sj- We will 
call a specific choice of contours weighted by these values P%, again leaving the dependence 
on u and v implicit. The choice u = \ and v = 1 again gives a particularly simple solution. 

Before turning to the extraction procedure, we may observe that the four-dimensional 
heptacuts do not suffice to extract information about the coefficients beyond O(e ). The 
problem is that we can find non-vanishing linear combinations of tensor integrals whose 
heptacut integrand vanishes identically for all six solutions. As a result, not only do integrals 
over all contours T^ a vanish, but even integrals constructed by multiplying the integrand 
by an arbitrary function of the remaining degree of freedom z would vanish. We call such 
linear combinations magic. Examples of magic combinations include, 

Mt = P%[1, 1] + S fP^ 2 [0, 1] - I (^) 2 P 2 *; 2 [0, 0] , 
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M 2 = P**[2, 2] + ^P 2 ** 2 [2, 1] + ^P 2 *; 2 [l, 2] - x (^) 2 P 2 r 2 [l, 1] , (4.26) 
where the abbreviated notation P 22 [m, n] is defined by, 

p** 2 K n ] = p 2 **[(^ • h) m (e 2 ■ hr\ . (4.27) 

The magic combinations do not vanish, but both coefficients of master integrals are of 0(e) 
after use of IBP reduction equations. 

V. INTEGRAL COEFFICIENTS 

With solutions to the constraint equations that also isolate specific basis integrals in hand, 
we can write down a procedure for computing the coefficients of the integrals in the master 
equation (14.11) . To do so, we apply the augmented heptacuts to the left-hand side of the 
master equation. The basic heptacut will break apart the two-loop amplitude into a product 
of six on-shell tree amplitudes, one for each vertex in the double box. We will be left with 
the integral over the z contour. On the right-hand side, we have the two basis integrals (14. 6 p 
chosen earlier. Here, apply the augmented heptacut, and perform all integrations. This 
gives us the relation, 

128 VW UJv^iz + xY } \} x 3 K) l6xs 3 12 f^ {1 ' J 2 ' J) 32 S ? 2 ^ 3 lj 

(5.1) 

In this equation, the product of amplitudes arises from a factor of a tree-level amplitude at 
each vertex of the double box with all seven propagators cut. 

As explained in the previous section, through a judicious choice of contours, we can make 
the coefficient of c 2 in this equation vanish, or alternatively the coefficient of c\ vanish. This 
would then allow us to solve for c\ and c 2 , respectively. We gave such choices in eqs. (I4.24p 
and (I4.25p . Using them, we can write an expression for ci, 





8u J Pl z[z + x) fJi J 


(5.2) 


and for c 2 , 




4si 2 u Jp 2 z{z + x) fJi 3 


(5.3) 
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The right-hand sides of these equations must be summed over possible helicity and particle- 
species assignments. The explicit integration is understood to be over the z component of P\ 
and P 2 respectively, with the integrations over the other a, and implicit in the solutions 
<Sj, and with the dependence of Pj on the parameters u and v left implicit. The formulae 
( 15. 2p and (15.31) represent the central result of this paper. They are valid for any gauge theory, 
and indeed for any amplitude satisfying the power-counting rules of gauge theory. With the 
notation, 

v\ Sl+S2 - S:i = v\ Si + v\ S2 - v\ 

we can write out these formulae more explicitly, 



is 3 ' 



(5.4) 



ci 



(v — 2u)i 
8^ 



Res - Y\Af ee (z) 

2 =0 z - LJ - 3 



IV 



Res 



S 1 +S 2 +S 3 +S i -S s -S 6 



m z=-x z + x 



I[Af-{z) 



3=1 



-- Res 



c 2 



z=-x z + x 



2u)i 



+ 



IX 



Res-ffAf c (z) 



Res 

4si2X 2 =-x z + X 



S1+S2+S3+S4 



4(1 + X) 
(6u — v)i 



<Ss +<Se 



Res-TTAf ee (z) 

4s 12 ux 2=0 2: ^ 



tree / 



+ 



Res 



5i+52+<S 3 +<S 4 



4si 2 mx 2 =-x z + x 
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tree / 



3i 



2s 12 (l + x) 



Res TTAf ee (2) 



i=i 



(5.5) 



These formulae are not manifestly independent of the choice of contour, but the constraint 
equations ensure that they are. We will see explicit examples in the next section. Of course, 
the independence of the final result of the choice of contour does not mean that the results at 
intermediate steps are independent; certain choices of contour may in fact simplify analytic 
or numerical calculations. We have already seen hints of this in the choices of P\ and P2, 
where some values of u and v will require evaluation of fewer contours, and hence possibly 
fewer numerical evaluations if the formulae (15. 2 p and (I5.3P are used in a numerical setting. 

At one loop, one can choose a basis so that integral coefficients are independent of the 
dimensional regulator e, and four-dimensional cuts suffice to compute all of them. (Com- 
puting the rational terms requires use of D-dimensional cuts.) At two loops, the coefficients 
of integral reductions, and hence generally of integrals in eq. ( I4.ip . will depend explicitly 
on e. In particular, C\ and c 2 above will depend explicitly on e. In general, this depen- 
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FIG. 5: Graphical representation of n*=i ^j ee {z) for the all-plus QCD amplitude 



dence cannot be extracted from four- dimensional heptacuts alone, because of the vanishing 
of magic combinations discussed in Section |TVJ We can also see the need for cuts beyond 
four dimensions, or considerably relaxing some of the heptacut conditions, by considering 

f +), computed in ref. 



the two-loop all-plus amplitude, A\ loop (+ 



531 ] . In this case, the 



product of tree amplitudes in eqs. (I5.2p and (15. 3p will necessarily vanish in four dimensions, 
because there is no assignment of internal helicities in fig. [5] that will leave all three-point 
amplitudes non-vanishing. The same observation still holds if we relax some of the cut 
conditions, examining hexacuts or pentacuts. 



VI. EXAMPLES 



In this section, we apply the formalism developed in previous sections to several examples 
of two-loop four-point amplitudes. We use the master formulae (15.21) and (15.31) to compute the 
coefficients to O(e ) of the two double box basis integrals, P 2 *2lXI an d ^jM^i'M- We consider 
three different contributions to four-gluon amplitudes in supersymmetric theories with M = 
4, 2, 1 supersymmetries: the s- and t-channel contributions to .A4~ loop (l~, 2", 3 + , 4 + ), and the 
s-channel contributions to A%~ loop (l~ , 2 + , 3~, 4 + ). (The t-channel contributions to the latter 
amplitudes can be obtained by relabeling the arguments of the s-channel contribution.) 

We will express the results as multiples of the tree-level four-point amplitudes, 

i(12) 3 



(23) (34) (41) ' 



(6.1) 
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and 



(6.2) 



2" 3+ 




1- 4+ 



FIG. 6: The labeling used in Section IVH here shown for the s-channel contribution to 
AMoqP(i_ 2- 3+,4+). 

In this section, it will be convenient to have a label for each cut propagator in the double 
box. Accordingly, we adopt a different labeling from previous sections. It is displayed in 

fig. El 



A. The s-channel contribution to A^~ loop (l ,2 ,3 + ,4 + ) 



For this contribution, shown above in fig. [61 the helicities of the external states allow only 
gluons to propagate in either loop. For this reason, we will get the same result independent of 
the number of supersymmetries. We find that for all six solutions to the on-shell equations, 

6 

]jAf ee = -ts 2 12 s 23 A^ ++ . (6.3) 

5=1 

We can then use eq. ( 15. 2 1) (or equivalently the first equation in eq. ( 15. 51) ) to obtain, 

= -sl 2 s 23 A^ ++ ; (6.4) 
and eq. ( 15. 3 p (or equivalently the second equation in eq. ( 15. 5p ) to obtain, 



s l2 ux 2s 12 ux suX 2s 12 ux, 
0. (6.5) 
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We see that the dependence on the parameters u and v has disappeared, as expected. 
In the M = 4 theory, these turn out to be the exact coefficients; in theories with fewer 
supersymmetries, there are additional terms of 0(e) in these coefficients. 



B. The t-channel contribution to A 



2- loop , 



1-2-3+,^ 




FIG. 7: The heptacut for the ^-channel contribution to A\ loop (l ,2 ,3 + ,4 + ). 

We turn next to the computation of the coefficients in the t-channel contribution to the 
same amplitude considered in the previous section. The heptacut for this contribution is 
shown in fig. [7J In applying the formulae for the coefficients, we have cyclicly permuted the 
external momentum arguments, (1,2,3,4) — > (4, 1,2,3), so that we must replace % ~ X 1 - 
Otherwise, they are of course unchanged. 

In this contribution, computing the required products of tree amplitudes is more involved, 
and the computation also requires sums over supermultiplets of states propagating in the 
loops. As an example, we work through the computation of the product in solution S 2 ■ 
We have two possible helicity assignments for the internal lines, shown in fig. [SJ For gluon 
internal lines, we multiply the amplitudes at the six vertices to obtain, 

A A for configuration A , 

(6.6) 

B 4 for configuration B , 
where 

A 4 = ([p4](lp 2 ) (q 2 l) [lq] (q qi ) [ qi q 2 ])\ 



6 1 

nttree, gluon _ x 

A i " "A 

7=1 
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(a) 



(b) 



FIG. 8: The two distinct assignments of internal helicities for the t-channel double-box contribu- 
tions to A^~ loop (l~, 2~, 3 + , 4 + ): (a) configuration A and (b) configuration B. 

5 4 = ( [4 Pi] (pi 1) (pa ft) [qp] (q ft) [ft ft] ) 4 , (6.7) 
A = [p4] [4pi] [pip] (pi 1) (lpa) (papi) (p 2 ft) (ft {Ipz) [lq] [qp] [pi] 
x(ftft) (ft 3) (3q) [ qi q 2 ] [q 2 2] [2q x ] , 

and the minus sign in eq. ( 16. 6 p comes from the factor of i in each y 4* ree ' § luor \ 

The helicity assignments of the internal lines allow only gluons to propagate in the right 
(q) loop, whereas the entire supersymmetric multiplet of states can propagate in the left (p) 
loop. For Af = 4 super Yang-Mills, the sum over states yields, 

(A + B) 4 



tree 
3 



M=4 j = l 

multinlet 



On the other hand, from refs. 



36j and {fill we know that in the j\f = 4 theory, 



A 



(6.8) 



s 2 



£ IK 

JV=4 j = l 
multiplet 



-is 12 s* A 1 ™ 



(6.9) 



S2 



As a calculational shortcut, we use the equality of the expressions in eqs. (16. 8 j) and (16. 9p 
to fix the relative sign of A and 5 in eq. ( 16. 7p . (Of course, the relative signs can also be 
determined a priori, without reference to results in the literature, by carefully tracking the 
direction — incoming or outgoing — of the momenta at a given vertex and using the analytic 
continuation rule that changing the sig n of a momentum, pi — > —p i: is effected by changing 
the sign of the holomorphic spinor [76j: Af — > — Af while Af — > Af.) One finds, 

A= [p4] (lp 2 ) (q 2 l) [lq] (q qi ) [ qi q 2 ] , 
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B = - [4p x ] (pi 1) (p 2 g 2 ) [gp] (ggi) [giga] • (6.10) 

Ref. [77j teaches us that the sum over the M = 4, 2, 1, multiplet of states is related to the 
M = 4 state sum via 

e rW = e rW. (o-u) 

SUSY j — 1 V / Af=4 j — 1 

multiplet multiplet 

so that the sum over the supersymmetric multiplet of states can be calculated from the 
gluonic contributions alone (indeed, recall that A and B in eq. (16.71) were obtained from the 
product of purely gluonic amplitudes corresponding to configurations A and B, respectively). 

We can simplify the expression for the ratio between the supersymmetric state sums in 
eq. (16. lip by factoring out as many common factors of A and B as possible (exploiting 
momentum conservation fully). Setting A = aF and B = (3F, for M = 4, 2, 1 the ratio 
appearing in eq. (16. lip simplifies to 

R = {^TW ( ~ 1Na) = (a + ( } 

where the last equality holds only for M = 4,2,1; it can be obtained by expanding the 



-a. 



numerator (a 4 + /3 4 ) (l — in eq. (I6.12p in (5 around 

In the case at hand, we can use momentum conservation (I — p 2 + g 2 and Pi = p — £4) to 
rewrite A and B as follows, 

A= [p4] (lp 2 ) (?2P2) [P2q] (qqi) [?i 92] , 

B = - [4p] (pi) (p 2 g 2 ) [gp] (ggi) [gife] , (6-14) 



and identify, 



a = (lp 2 ) [p 2 g] , 

P = -[?p](pl> , (6.15) 
F = [p4] (g 2 p 2 ) (ggi) [gi g 2 ] ■ 



Momentum conservation implies that a + /9 = — (14) [4g], and thus, 

Q = (lpi) [pi g] = (lpi) [pi 3] 
a + (14) [4 g] (14) [43] ' 



(6.16) 
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where the second equality uses the proportionality of antiholomorphic spinors, X q oc A3. 
(This proportionality holds only for some of the other six solutions <Sj in addition to 1S2.) 
The ratio thus simplifies to, 

*- i+ *MWf) +l *MWFf- (6 - 17) 

We can solve for the explicit values of the cut momenta using the parametrization ( 13 .4p with 
the external momenta cyclicly permuted (for the t channel configuration), and using the 
on-shell values defining <S 2 given in eqs. ( 13.7ft and (13.101) . We find, 

py=p"-Aff= , SuZ r . (4~ 1 7 M |1~), (6.18) 
2 (4 3) [3 1] V ' 

so that, 

(1ft) M = (1-1 Y 13") iv - 2 (4 S 3^3 1} (11 r ]3-> (41 7m 11-) 

(41) , . 

= -py s ^, (6.19) 

and thus, 

<^M = X *. (6-20) 

(1 4) [4 3] A K ' 

This gives us our final form for the ratio, 



R= 1 + (4 - Af) X z + {4 - Af) X z 2 , (6.21) 



and for the product of tree amplitudes, 

6 



E U A T =-is 12 sl 3 A^ + Jl + (4-Af) X z + (A-N-)x 2 ^). (6.22) 



SUSY j — 1 
multiplet 



In this solution to the heptacut equations, the supersymmetric multiplet runs only in 
one of the loops. In other solutions (in particular, Sq), the multiplet can run in both loops. 



The treatment of this case is similar but more elaborate. It turns out [77j] that the sum 
over the multiplet can again be evaluated purely from the gluonic contributions. The main 
difference is that in this case there are three gluonic contributions A 4 , B A , C A (compared to 
the two in eq. (16.71) ). One can again fix the relative sign of B and C by insisting that the 
N = 4 supersymmetric result — ( A+B + c ^ be equal to eq. (16 .9p . and from the obvious analog 
of eq. ( 16.1 ip one then finds the results for the supermultiplet sums for Af = 4, 2, 1, 0. These 
expressions can again be simplified as above. 
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Summing over all six solutions, and plugging the result into our master formula? (I5.2p 
and f)5.3p . taken with u = | and v — 1, we find 



Cl = -Sl2« 

3 



23^--++ (l + ^(1 - <M(4 - AT)! X ( X + 1)^ 



c 2 



2 /i tree 



valid for M = 4, 2, 1. 



(6.23) 



C. The s-channel contribution to A^~ loop (l ,2 + ,3 ,4 



)+ 




l + 



FIG. 9: The heptacut for the s-channel contribution to A\ loop (l ,2 + ,3 ,4 + ). 

The heptacut for the s-channel contribution to 7i4~ loop (l~, 2 + , 3~, 4 + ) is shown in fig. EO 
We will evaluate this contribution in two different ways, illustrating both the result's inde- 
pendence of the precise choice of contour, and also illustrating the potential advantages of 
a judicious choice of contour in a given calculation. 

Rather than using our master formula? (I5.2p and (15. 3p . let us evaluate the augmented 
heptacut integral for a general contour, before imposing the constraint equations. Adding 
up the contributions from all six solutions, we find 

dz 



6 



1 Jvi A* + x) fJi 

-isi 2 s 23^4-+-+ 



E 



X 



(4 -AO 



Ol,6 — a 3,6 ~~ a 2,5 



(x 



+ (1 - §M 



X 



4-Af 



(x + i) 



4-JV 



- 1 



Qj,3 ~ a 2,3 a l,4 — a 2,4 



X 



.(6.24) 
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In this expression, we need to impose the constraint equations in order to restrict the 
evaluation to a valid contour; and then we would seek to project onto each basis integral 
in turn. Now, suppose we can find a pair of solutions to the constraint equations which 
projects onto the first or second basis integral, respectively, and in addition, satisfies 0^3 — 
^2,3 + Qi,4 — o 2)4 = 0. Using such a contour would set the second line of eq. (16.241) equal to 
zero and therefore produce a particularly simple algebraic expression for c\ and c 2 directly, 
without need for additional simplification. Choosing u = | and v — 1 in Pi and P 2 gives 
such a contour. 

This gives us the results, 

02 = -ls 12 s^_ + ±^, (6.25) 

valid for M = 4, 2, 1. The t-channel contribution can be obtained by exchanging s 12 -H- s 2 3 
and x -> X -1 - 

If we compare the expressions obtained above for the coefficients q to those obtained 
using the choice suggested in Section [V] u—\ and v — 1, we find that the expressions are 
equal by virtue of the identity 

valid for M = 4, 2, 1. This identity can of course easily be proven without reference to the 
current discussion, but the point we wish to emphasize is that the flexibility in choosing 
contours suggests certain algebraic simplifications which are not immediately obvious. 

The double box coefficients given in eqs. (16 A\ 16.51 16.231 16. 25ft agree with the O(e ) terms 
of the corresponding coefficients, supplied to us by Lance Dixon 78], in the amplitudes 
computed by Bern, De Freitas, and Dixon |s4]. 



VII. CONCLUSIONS 



In this paper, we have taken the first step to extending the maximal generalized uni- 
tarity method to two loops. Cutting propagators can be viewed as deforming the original 
real loop-momentum contours of integration to contours encircling the global poles of the 
integrand. At two loops, there is a variety of such poles. We can evaluate the integral along 
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many different linear combinations of these contours. However, our choices are restricted by 
the requirement that the evaluation along any contour respect the vanishing of certain inser- 
tions of Levi-Civita symbols, as well as of total derivatives arising from integration-by-parts 
identities. We derived the corresponding constraint equations for the massless double box, 
and showed how to use their solutions to obtain simple formulae, eqs. (15. 2p and ( 15. 3p . for the 
coefficients of the two double box basis integrals to leading order in the dimensional regu- 
lator e. To derive these equations, we adopted a parametrization of the loop momenta and 
solved explicitly for the maximal cuts, a heptacut in our case, and identified the additional 
poles present in the remaining degree of freedom. 

We expect that the approach given in this paper — parametrize the basis integrals; solve 
the on-shell equations; identify the poles in the remaining degrees of freedom; impose all 
constraint equations — will apply to the full set of integrals required for two-loop amplitudes, 
both to the four-dimensional cuts considered here, and more generally to the .D- dimensional 
cuts required for a complete calculation of the amplitude. 
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